Understanding the different scaling behavior in various shell models proposed for 

turbulent thermal convection 
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Different scaling behavior has been reported in various shell models proposed for turbulent thermal 
convection. In this paper, we show that buoyancy is not always relevant to the statistical properties 
of these shell models even though there is an explicit coupling between velocity and temperature 
in the equations of motion. When buoyancy is relevant (irrelevant) to the statistical properties, 
the scaling behavior is Bolgiano-Obukhov (Kolmogorov) plus intermittency corrections. We show 
that the intermittency corrections of temperature could be solely attributed to fluctuations in the 
entropy transfer rate when buoyancy is relevant but due to fluctuations in both energy and entropy 
transfer rates when buoyancy is irrelevant. This difference can be used as a criterion to distinguish 
whether temperature is behaving as an active or a passive scalar. 

PACS numbers: 47.27-i,47.27te 



I. INTRODUCTION 



Turbulent thermal convection is a_problem of great re- 
search interest (see, for example, for a review). One 
interesting issue is to understand the scaling behavior 
of the velocity and temperature fluctuations. Turbulent 
thermal convection is often investigated experimentally 
in Rayleigh-Benard convection cells, which are closed 
cells of fluid heated from below and cooled on the top. 
Such confined turbulent convective flows are highly inho- 
mogeneous as thermal and viscous boundary layers are 
formed near the top and bottom of the cell. Scaling laws 
for the central region of such confined turbulent ther- 
mal convection have been put forth and shown to be in 
good agreement with the existing experimental measure- 
ments [3]. On the other hand, shell models focussing on 
the energy cascade process have been studied intensively 
and proved to be useful for understanding the scaling 
behavior of velocity fluctuations in inertia-driven turbu- 
lence (see, for example, 4] for a review). It is thus natural 
to also construct shell models for turbulent thermal con- 
vection. Shell models are, by construction, boundary-free 
and thus shell models for turbulent thermal convection 
are necessarily models of homogeneous turbulent thermal 
convection. It is known that the presence of boundaries 
generates coherent structures such as plumes and a large- 
scale mean flow in confined turbulent thermal convection, 
and these coherent structures can affect the scaling be- 
havior [3|. Thus, scaling behavior in confined turbulent 
thermal convection and scaling behavior in homogeneous 
turbulent thermal convection as studied in shell models 
can be different. 

Several shell models for turbulent thermal convec- 
tion have been proposed and different scaling behav- 
ior reported. Specifically, Bolgiano-Obukhov (BO) scal- 
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ing Q plus intermittency corrections has been reported 
in the shell model constructed by Brandenburg [6| and 
also in the modified model by Suzuki and Toh Q 
for some parameter range. On the other hand, Kol- 
mogorov 1941 (K41) scaling ^] plus intermittency cor- 
rections has been reported by Jiang and Liu Q using a 
shell model extended from the Gledzer-Ohkitani-Yamada 
(GOY) model which we shall denote as the GOYT 
model. In this paper, we show that buoyancy is not al- 
ways significant and directly relevant to the statistical 
properties even though there is an explicit coupling be- 
tween velocity and temperature in the equations of mo- 
tion in all these shell models. We clarify that the two 
different types of scaling behavior reported correspond 
respectively to the case when buoyancy is relevant to the 
statistical properties and the case when it is not. Specif- 
ically, the scaling behavior is BO plus intermittency cor- 
rections when buoyancy is relevant, and K41 plus inter- 
mittency corrections (as one would expect for tempera- 
ture behaving as a passive scalar) when buoyancy is ir- 
relevant. We show that the intermittency corrections of 
temperature could be solely attributed to fluctuations in 
the entropy transfer rate when buoyancy is relevant but 
due to fluctuations in both energy and entropy transfer 
rates when buoyancy is irrelevant. This difference might 
be used as a criterion to distinguish whether temperature 
is behaving as an active or a passive scalar. 



II. SHELL MODELS PROPOSED FOR 
TURBULENT THERMAL CONVECTION 

Two classes of shell models have been proposed for 
studying turbulent thermal convection. The first class 
consists of the shell model proposed by Brandenburg ^ 
and its modified versions The other class consists 
of the GOYT model, the shell model extended from the 
GOY model [OJ and the SabraT model [llj from the Sabra 
model The Sabra model [l2[ was proposed to elim- 
inate some undesirable periodic oscillations in the GOY 
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model, and have essentially the same scaling behavior as 
the GOY model. The scaling behavior in the first class 
of shell models is BO plus corrections in some parameter 
range while the scaling behavior in the second class of 
shell models is always K41 plus corrections. In this pa- 
per, we focus on two shell models, one from each class. 
The first one, denoted as the Brandenburg model, is the 
modified model proposed by Suzuki and Toh [7,] without 
the drag term. The second is the SabraT model. 

The basic idea of a shell model is to consider variables 
in discrete "shells" in Fourier /c-space, and construct a 
set of ordinary differential equations for these variables 
per shell. For shell models for turbulent thermal convec- 
tion, there are two variables, the velocity and tempera- 
ture variables, u„ and 0„. They can be roughly thought 
of as the Fourier transforms of the velocity and tempera- 
ture fields with wavevector fc, whose magnitude satisfies 
kn l£ \k\ < kn+i- Here, fc„ = 2"fco is the wavenumber 
of the nth shell, with < n < — 1, and fco = 1 is 
the wavenumber corresponding to the largest scale in the 
system. The equations of motion for u„ and 0„ are: 



^ - I (k 
dt ~ " 

do 

= l0{kn) - Kk^dn + fn 



ag6„ 



(1) 
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where /„ is the forcing term acting only on the first few 
shells. The nonlinear terms luikn) and Ie{kn) are taken 
to couple quadratically with the nearest shells and some- 
times also the next nearest shells, and are constructed 
to satisfy two conservation laws of energy and entropy 
(proportional to |0n|") in the limit of i/ ^ and k ^ 0: 
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As a result, the nonlinear terms u^Iuikn) and 9'^Ig{kn) 
should have a fluxlike form such that the evolution equa- 
tions of energy and entropy in the nth shell are: 



d_ 

dt 

d_ 

di 



Fuikn) - Fu{k„+i) - iykl\u„\^ (5) 
+agRe{unOn} 



The fluxes Fu{kn) and Fg{kn) are respectively the rates 
of energy and entropy transfer from the n — 1th shell to 
the nth shell. 

In the Brandenburg model, w„ and 9n are real variables 
with @,0]: 



(kn) 

F^ikn) 

FFikn) 



dkn{Un-lOn^l - 2w„6'„+i) 

+bk 

n ('^n^n — 1 '^'^n-\-l^ n-\-l ) (8) 

[aun-i + hun)k (9) 
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where a, &, a and b are positive parameters. In the 
SabraT model, m„ and On are complex variables with [ll[: 

luikn) = iknMUn+iUn+2 ~ -jKi^-^Un+l 

1-5 

H —Un-lUn^2), (H) 

Igikn) = ikn{aiK^+idn+2 + Ol2Un+20n+l 

- P2Un+lO'^^i 
-JlUn-lOn-2 - l2Un~26n-l) (12) 

Fu{kn) = AIm[fc„_iu*_iU*u„+i 

+ {l - ^)fc„_2<-2<-l^n] (13) 
Fg(kn) = lTa[^l{knUn-lB.n-20*n + kn+lU.nOn^lO*n+i) 

-f32knU*n+ld,i^l0n + 12k ](14) 

The parameters ai,2, /3i,2 and 71^2 are determined by 
ai = 4t , /3i = 1 — (5 — 2t , 71 = — r, 
a2 = 2 - 4t , 02 = 1 - 2t , 72 = r - (15) 

with three free parameters A, S and r. In particular, we 
fix A = 2 and r = 0.7 and vary S. The value (5 = 1 is the 
boundary value separating two families of Sabra model: a 
family of three-dimensional-like models for < 5 < 1 and 
a family of two-dimensional- like models for 1 < (5 < 2. 
We focus on < 5 < 1 in this paper. 

We study the scaling behavior of the velocity and tem- 
perature structure functions, (|wn|^) and (|6'„|^), with 
scaling exponents Cp and defined by: 



iKn^.k-i" 



(Kn-k-i^ 



(16) 



(7) 



where (. . .} denotes a time average. The K41 scaling 
would be characterized hy (p = £,p = p/3 while the BO 
scaling by (^p = Bp/ 5 and £,p = p/5. In our numerical 
calculations, we integrate the equations of motion using 
fourth order Runge Kutta method with an initial condi- 
tion oi Un = On — except for a small perturbation of 
9n at intermediate values of n. The Brandenburg model 
is forced with /„ = f5n,o where / is a uniform random 
noise while the SabraT model is forced with a Gaussian 
time-correlated noise acting on n = 3 and 4 only p^ . 
For the results presented in this work, we summarize the 
parameters used in Table HI 

In the Brandenburg model, the scaling behavior de- 
pends on the relative magnitudes of the parameters a 
and b, as reported in earher studies @. When 6/a is 
larger than some critical value of about 2, the scaling 
exponents and ^p are given by the BO values plus cor- 
rections. The scaling behavior improves with b/a. On 
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TABLE I; Values of the parameters used for the results pre- 
sented. 



Brandenburg model 


a 


b 


a and b 


V 






iV 


0.01 


1 


1 


5 X 10"" 


5 X 10"^^ 


1 


32 


0.31 


0.6 


1 


5 X lO"'' 


5 X lO"'' 


1 


25 


SabraT model 


S 


A 


r 


V 


K 


ag 


iV 


0.5 


2 


0.7 


10"* 


10"* 


1 


23 


0.8 


2 


0.7 


10"* 


10"* 


1 


23 



the other hand, when hja is smaller but close to the crit- 
ical value, the scaling exponents C^p and are the same 
as those obtained in the case of passive scalar advection 
in which the couphng term agBn with temperature in 
the velocity equation of motion is replaced by a random 
forcing at n — 0. This indicates that buoyancy does 
not play a part in the statistical properties in this case. 
The scaling exponents for hja = 100 and h/a = 1.94 are 
shown respectively in Figs. [1] and [D For even smaller 
values of 6/a, further away from the critical value, the 
system is not chaotic, and in most of the shells the solu- 
tion is given instantaneously by the fixed-point solution 
of Un — Akn^^^ and 0„ = Bkn^^^ , which holds exactly 
in the limit of large N and v = k. ~ ag = Q. 




FIG. 1: The scaling exponents C,p (squares) and (circles) 
for Brandenburg model with a — 0.01 and b — 1. The error 
increases with p and the largest errors are shown. Comparing 
with the two solid lines of slopes 1/5 and 3/5 shown, it can 
be seen that the scaling behavior is BO with corrections. 



For the SabraT model, we find that the values of Cp 
remain the same as those in the Sabra model without 
the coupling term agOn for all the values of 6 studied, 
again indicating that buoyancy does not play a role in de- 
termining the statistical properties in the SabraT model 
for < (5 < 1. The precise values of (p depend on S, 
as was reported in the GOY model [Ts'l . In Fig. [H we 
present the results for and for S = 0.5, a con- 




FIG. 2: Same as Fig. [T] for a = 0.31 and b = 0.6. The solid 
line shown has slope 1/3. 

ventional value at which the model conserves helicity 
in the inviscid limit [l^. In this case, the values of 
(p are well described by the She-Leveque result [3] of 
= p/9 + 2[1 - (2/3)P/^], as was also reported ^ for 
the GOYT model with S = 0.5. 




FIG. 3: Same as Fig. [2] for the SabraT model with 6 = 0.5. 
The solid line shown has a slope of 1/3 while the dashed line 
is the She-Leveque result [13]. 



III. THE BUOYANCY SCALE 

In this section, we discuss how to determine whether 
buoyancy is relevant or not in determining the statistical 
properties. Consider Eq. which is the energy budget. 
The third term of the right hand side is the rate of energy 
dissipation in the nth shell due to viscosity while the 
last term is the power injected into the nth shell by the 
buoyancy forces. It is thus reasonable to take buoyancy 
to be significant in the nth shell if 

\ag{Re{un0:})\>e (17) 

where e = ^^X^n ^n(l^nP) the average energy dissipa- 
tion rate. We denote the scale at which the equality sign 
in Eq. ([T7]) holds to be the buoyancy scale fc„. . Hence 
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buoyancy is relevant and significant for n < n* and ir- 
relevant or insignificant for n > n*. It is easy to show 
that for Un and 9n satisfying exactly K41 or BO scaling, 
fc„. = 1/Lb, where Lb = e'/^x~^'''{ag)~^''^ is the Bol- 
giano length [isj and x is the average thermal or entropy 
dissipation rate given by x = kX)™ ^n(l^nl^)- 

As shown in Figs. [3] and [51 we find that Eq. ([TT]) is 
satisfied for most of the shells only in the Brandenburg 
model with h/a larger than the critical value. When h/a 
is smaller than the critical value, buoyancy is insignificant 
in all except the largest shells. For the SabraT model, 
we find that buoyancy is insignificant in all except the 
largest shells for all the values of 5 studied. The results 
for 5 — 0.5 and 5 — 0.8 are shown in Fig. [6l 
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FIG. 4: Comparison of \ag{un9n)\ (circles) with e (solid line) 
in each shell for the Brandenburg model with a large value of 
h/a = 100. 



The two different scaling behavior can be understood by 
studying the evolution equations of energy and entropy. 
In the intermediate range where external forcing is not 
acting and where energy and entropy dissipation rates 
are both small, Eqs. (O and ([6]) can be approximately 
written as: 



Fu{K) 



Fu{k„+i) + agRe{un9*J « (18) 
Feikn) - FeiK+i) « (19) 



From Eq. (|19p . Fg{kn) is independent of /c„ in the in- 
termediate range, implying that there is an entropy 
cascade. From Eq. ^5)) . we see that agRe{unO*^} is 
comparable with Fu when buoyancy is significant, and 
Fu{kn) — Fu{kn+i) « when buoyancy is insignificant. 
Thus when buoyancy is insignificant, there is also an en- 
ergy cascade as in the usual inertia-driven turbulence. 
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FIG. 6; Comparison of \ag{Re{un9^})\ (circles) with e (solid 
line) in each shell for the SabraT model with S — 0.5 in the 
top panel and 5 = 0.8 in the bottom panel. 



FIG. 5: Same as Fig. |4]for a small value of b/a ~ 1.9. 

One naturally expects different scaling behavior when 
buoyancy is significant and when it is not. In this sense, 
it is not puzzling that different scaling behavior was re- 
ported in the the various shell models proposed. Indeed 
we find BO scaling plus corrections when buoyancy is sig- 
nificant and K41 scaling plus correction when it is not. 



In the case when buoyancy is significant, there is only 
the cascade of entropy. As a result, one expects the statis- 
tical properties to be controlled by the entropy cascade. 
Specifically, one expects fl^ the statistical properties of 
Un and On to be determined solely by Fg, ag, and fc„: 

\un\ = 0„(ag)2/5|^^,(fc„)|i/5fc-3/5 (20) 
l^nl = </>e(aff)-^/5|Fe(fc„)|2/5fc-i/5 (21) 
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where and (fig are dimensionless random variables that 
are independent of and statistically independent of 
Fg{kn). On the other hand, when buoyancy is insignifi- 
cant, there is also the energy cascade. Thus one expects 
the statistical properties of u„ and On to be determined 
by Fu, Fg and fc„: 

\Un\ = V«|J^u(fc„)r/'fc-l/3 (22) 

\0n\ = V'e|J^.(A:„)ri/^|Fe(fc„)|i/2fc-i/3 (23) 

where "011 and ipg are dimensionless random variables that 
are independent of kn and statistically independent of 
Fu(kn) and Fg{kn)- Hence we have 

(knD - (|F,(fc„)r/^)A:-3p/5 (24) 
{\0nn - (|F,(fc„)pP/5)fc-p/5 (25) 

when buoyancy is significant and 

- {\FMf^')K'^' (26) 
(Kn ^ {\F^{kn)rP^''\Fe{knW/^)k-P/' (27) 

when it is not. Equations and (PSjl . and Eqs. ([^ 
and (P7|) thus respectively give BO and K41 scaling plus 
intermittency corrections for the case when buoyancy is 
significant and when it is not, just as what was found 
numerically. Moreover, when buoyancy is significant, the 
intermittency corrections are solely due to fluctuations 
in Fg while in the case when buoyancy is insignificant, 
the intermittency corrections are due to fluctuations in 
both and Fu and Fg. We have checked and verified [l^ 
Eqs. p4| and ((25l) . 




FIG. 7: The scaling exponents ^* of tlie conditional tempera- 
ture structure functions (|Sn|^ | Fe — x) for the Brandenburg 
model with small b/a (squares) and large h/a (circles), and the 
SabraT model with 5 — 0.5 (triangles). The error increases 
with p and the largest errors are shown. Two solid lines with 
slopes 1/3 and 1/5 are shown. 

Our work shows that the mere presence of a coupling 
term between velocity and temperature in the equations 
of motion does not automatically imply that buoyancy 
is significant and affects the statistical properties. This 



leads to the question: How can one tell whether tempera- 
ture is behaving as an active or a passive scalar in models 
for turbulent thermal convection? For shell models, one 
can use Eq. (|17p . If Eq. ([TT]) is satisfied in most shells 
then buoyancy is significant and temperature is active 
otherwise temperature would behave as a passive scalar, 
ft would also be useful to have some other criterion that 
involves directly the statistical features of temperature. 
Equations (pO| and (|2T|) imply that when buoyancy is sig- 
nificant, the conditional statistics of u„ and 9n at fixed 
values of Fg would have simple BO scaling with no cor- 
rections [l6| . On the other hand, this is not true when 
buoyancy is insignificant; instead Eqs. ((22)) and ([23]) indi- 
cate that the conditional statistics of m„ and 0„ at fixed 
values of Fg continue to deviate from simple K41 scaling. 
Hence one can study the conditional statistics of temper- 
ature at fixed values of the entropy transfer rate. If these 
conditional statistics are described by simple scaling then 
temperature is behaving as an active scalar. Otherwise 
if the conditional statistics remain anomalous then tem- 
perature is behaving as a passive scalar. To check this 
idea, we calculate the the conditional temperature struc- 
ture functions at fixed values of entropy transfer rate and 
their scaling exponents 

iKl" \ Fg = x) k~^'- (28) 

for the SabraT model the Brandenburg model for both 
small and large values of b/a. The results of £,* do not 
depend on x and are shown in Fig. [Tj It can be seen that 
for Brandenburg model with large b/a, ^* are indeed well 
described by the BO values of p/5. Also, as expected, 
for both the Brandenburg model with small b/a and the 
SabraT model, ^*'s continue to deviate from the K41 
values of p/3. 



IV. CONCLUSIONS 

Various shell models have been proposed for turbu- 
lent thermal convection. K41 scaling plus corrections has 
been reported in most of these models while BO scaling 
plus intermittency corrections is reported in the Bran- 
denburg model with suitable parameters. In this paper, 
we have shown that buoyancy is not always significant 
and relevant to the statistical properties in these shell 
models even though there is an explicit coupling term 
with temperature in the equation of motion for velocity. 
We have further clarified that BO scaling plus correc- 
tions would be observed only in the shell models in which 
buoyancy is significant. For shell models in which buoy- 
ancy is insignificant, the statistical properties remain the 
same as in the case in which the coupling term with tem- 
perature is absent. We have argued that the statistics 
properties are controlled solely by the cascade of entropy 
when buoyancy is significant but controlled by both the 
cascades of energy and entropy when buoyancy is not 
significant, and shown how this leads to the two different 
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scaling behavior in the two cases. We have further shown 
that the intermittency corrections are solely attributed 
to fluctuations of the entropy transfer rate when buoy- 
ancy is significant but are caused by fluctuations of both 
the energy and entropy transfer rate when buoyancy is 
insignificant. As a result, the conditional temperature 
structure functions at fixed entropy transfer rate would 
have simple scaling when buoyancy is significant but re- 
main anomalous when buoyancy is insignificant. We have 
demonstrated how this feature can be used as a criterion 



to distinguish whether temperature is acting as an active 
or a passive scalar. 
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